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Primordial nucleosynthesis provides a probe of the Universe during its early evolution. 
Given the progress exploring the constituents, structure, and recent evolution of the 
Universe, it is timely to review the status of Big Bang Nucleosynthesis (BBN) and to 
confront its predictions, and the constraints which emerge from them, with those derived 
from independent observations of the Universe at much later epochs in its evolution. Fol- 
lowing an overview of the key physics controlling element synthesis in the early Universe, 
the predictions of the standard models of cosmology and particle physics are presented, 
along with those from some non-standard models. The observational data used to infer 
the primordial abundances are described, with an emphasis on the distinction between 
precision and accuracy. These relic abundances are compared with predictions, testing 
the internal consistency of BBN and enabling a comparison of the BBN constraints with 
those derived from the WMAP Cosmic Background Radiation data. Emerging from 
these comparisons is a successful standard model along with constraints on (or hints of) 
physics beyond the standard models of particle physics and of cosmology. 

Keywords: Big Bang Nucleosynthesis; Early Universe Expansion Rate; Neutrino Asym- 
metry; Abundances of D, ''He, ^He, ^Li. 

1. Introduction 

As the Universe evolved from its early, hot, dense beginnings (the "Big Bang") to 
its present, cold, dilute state, it passed through a brief epoch when the tempera- 
ture (average thermal energy) and density of its nucleon component were such that 
nuclear reactions building complex nuclei could occur. Because the nucleon content 
of the Universe is small (in a sense to be described below) and because the Uni- 
verse evolved through this epoch very rapidly, only the lightest nuclides (D, ^He, 
"^He, and ^Li) could be synthesized in astrophysically interesting abundances. The 
relic abundances of these nuclides provide probes of conditions and contents of the 
Universe at a very early epoch in its evolution (the first few minutes) otherwise 
hidden from our view. The standard model of Cosmology subsumes the standard 
model of particle physics {e.g., three families of very light, left-handed neutrinos 
along with their right-handed antineutrinos) and uses General Relativity {e.g., the 
Friedman equation) to track the time-evolution of the universal expansion rate and 
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its matter and radiation contents. While nuclear reactions among the nucleons are 
always occurring in the early Universe, Big Bang Nucleosynthesis (BBN) begins in 
earnest when the Universe is a few minutes old and it ends less than a half hour later 
when nuclear reactions are quenched by low temperatures and densities. The BBN 
abundances depend on the conditions (temperature, nucleon density, expansion 
rate, neutrino content and neutrino-antineutrino asymmetry, etc.) at those times 
and are largely independent of the detailed processes which established them. As 
a consequence, BBN can test and constrain the parameters of the standard model 
(SBBN), as well as probe any non-standard physics/cosmology which changes those 
conditions. 

The relic abundances of the light nuclides synthesized in BBN depend on the 
competition between the nucleon density-dependent nuclear reaction rates and the 
universal expansion rate. In addition, while all primordial abundances depend to 
some degree on the initial (when BBN begins) ratio of neutrons to protons, the ^He 
abundance is largely fixed by this ratio, which is determined by the competition 
between the weak interaction rates and the universal expansion rate, along with 
the magnitude of any Ue — De asymmetry.'* To summarize, in its simplest version 
BBN depends on three unknown parameters: the baryon asymmetry; the lepton 
asymmetry; the universal expansion rate. These parameters are quantified next. 



1.1. Baryon Asymmetry - Nucleon Abundance 

In the very early universe baryon-antibaryon pairs (quark-antiquark pairs) were as 
abundant as radiation {e.g., photons). As the Universe expanded and cooled, the 
pairs annihilated, leaving behind any baryon excess established during the earlier 
evolution of the Universe^. Subsequently, the number of baryons in a comoving vol- 
ume of the Universe is preserved. After pairs annihilate, when the temperature 
(in energy units) drops below the electron mass, the number of Cosmic Background 
Radiation (CBR) photons in a comoving volume is also preserved. As a result, it is 
useful (and conventional) to measure the universal baryon asymmetry by compar- 
ing the number of {excess) baryons to the number of photons in a comoving volume 
(post-e^ annihilation). This ratio defines the baryon abundance parameter tjb, 

,3 , (1) 

As will be seen from BBN, and as is confirmed by a variety of independent (non- 
BBN), astrophysical and cosmological data, tjb is very small. As a result, it is 
convenient to introduce jjio = 10^°?7b and to use it as one of the adjustable pa- 
rameters for BBN. An equivalent measure of the baryon density is provided by the 
baryon density parameter, JIb , the ratio (at present) of the baryon mass density to 



*A lepton asymmetry much larger than the baryon asymmetry (which is very small; see §1.1 below) 
would have to reside in the neutrinos since charge neutrality ensures that the electron-positron 
asymmetry is comparable to the baryon asymmetry. 
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the critical density. In terms of the present vahie of the Hubble parameter (see §1.2 
below), Hq = lOO/i km Mpc^^, these two measures are related by 

7710 = 10i°(nB/n^)o = 27417b/i^. (2) 

Note that the subscript refers to the present epoch (redshift z = 0). 

From a variety of non-BBN cosmological observations whose accuracy is dom- 
inated by the very precise CBR temperature fluctuation data from WMAP^, the 
baryon abundance parameter is limited to a narrow range centered near 7710 « 6. As 
a result, while the behavior of the BBN-predicted relic abundances will be described 
qualitatively as functions of tjb, for quantitative comparisons the results presented 
here will focus on the limited interval 4 < 7710 < 8. As will be seen below (§2.2), over 
this range there are very simple, yet accurate, analj^tic fits to the BBN-predicted 
primordial abundances. 



1.2. The Expansion Rate At BBN 

For the standard model of cosmology, the Friedman equation relates the expansion 
rate, quantified by the Hubble parameter (H), to the matter-radiation content of 
the Universe. 

= — GnPtot, (3) 

where Gn is Newton's gravitational constant. During the early evolution of the 

Universe the total density, ptoTi is dominated by "radiation" (i.e., by the con- 
tributions from massless and/or extremely relativistic particles). During radiation 
dominated epochs (RD), the age of the Universe (t) and the Hubble parameter are 
simply related by {Ht)RB = 1/2. 

Prior to BBN, at a temperature of a few MeV, the standard model of par- 
ticle physics determines that the relativistic particle content consists of photons, 

pairs and three flavors of left-handed (i.e., one helicity state) neutrinos (along 
with their right-handed, antineutrinos; N,^ = 3). With all chemical potentials set 
to zero (very small lepton asymmetry) the energy density of these constituents in 
thermal equilibrium is 

43 

PTOT = P7 + Pe + ^Pv = yP7, (4) 

where is the energy density in the CBR photons (which have redshifted to become 
the CBR photons observed today at a temperature of 2.7K). In this case (SBBN: 
Ni, = 3), the time-temperature relation derived from the Friedman equation is. 

Pre - e=^ annihilation : t = 0.738 MeV^ s. (5) 

In SBBN it is usually assumed that the neutrinos are fully decoupled prior to 
annihilation; if so, they don't share in the energy transferred from the annihi- 
lating pairs to the CBR photons. In this very good approximation, the pho- 
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tons are hotter than the neutrinos in the post-e^ annihilation universe by a factor 
T^/T^ = (11/4)1/3, and the total energy density is 

Ptot = Pf+ = 1.68p^, (6) 

corresponding to a modified time-temperature relation, 

Post - e='= annihilation : t = 1.32 MeV^ s. (7) 

Quite generally, new physics beyond the standard models of cosmology or parti- 
cle physics could lead to a non-standard, early Universe expansion rate {H'), whose 
ratio to the standard rate {H) may be parameterized by an expansion rate factor 
S, 

H^H' = SH. (8) 

A non-standard expansion rate might originate from modifications to the 3-1-1 di- 
mensional Friedman equation as in some higher dimensional models^, or from a 
change in the strength of gravity^. Different gravitational couplings for fermions 
and bosons^ would have similar effects. Alternatively, changing the particle popu- 
lation in early Universe will modify the energy density temperature relation, also 
leading, through cq. 3, to 5 7^ 1. While these different mechanisms for implement- 
ing a non-standard expansion rate are not necessarily equivalent, specific models 
generally lead to specific predictions for S. 

Consider, for example, the case of a non-standard energy density. 

Pk-^p'k = S^PR, (9) 

where = Pr + px and X identifies the non-standard component. With the 
restriction that the X are relativistic, this extra component, non-interacting at 
annihilation, behaves as would an additional neutrino flavor. It must be empha- 
sized that X is not restricted to additional flavors of active or sterile neutrinos. In 
this class of models S is constant prior to annihilation and it is convenient (and 
conventional) to account for the extra contribution to the standard-model energy 
density by normalizing it to that of an "equivalent" neutrino flavor^, so that 

7 

Px = AN^p^ = -AN^p^. (10) 
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For this case, 

S = Spre = (1 + ^^^-J • (11) 

In another class of non-standard models the early Universe is heated by the 
decay of a massive particle, produced earlier in the evolution"^. If the Universe is 
heated to a temperature which is too low to (re)populate a thermal spectrum of 
the standard neutrinos (Trh ^ 7 MeV), the effective number of neutrino flavors 
contributing to the total energy density is < 3, resulting in AA^,^ < and S < 1. 
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Since the expansion rate is more fundamental than is AiV,y, BBN for models with 
non-standard expansion rates will be parameterized using S (but, for comparison, 
the corresponding value of AN^, from eq. 11 will often be given for comparison). 
The simple, analytic fits to BBN presented below (§2.2) are quite accurate for 
0.85 <S< 1.15, corresponding to -1.7 < AN^ < 2.0 

1.3. Neutrino Asymmetry 

The baryon asymmetry of the Universe, quantified by r]B, is very small. If, as ex- 
pected in the ciirrcntly most popular particle physics models, the universal Icpton 
and baryon numbers are comparable, then any asymmetry between neutrinos and 
antineutrinos ("neutrino degeneracy") will be far too small to have a noticeable 
effect on BBN. However, it is possible that the baryon and lepton asymmetries 
are disconnected and that the Icpton (neutrino) asymmetry could be large enough 
to perturb the SBBN predictions. In analogy with rjB which quantifies the baryon 
asymmetry, the lepton (neutrino) asymmetry, L = = SaL^^ , may be quantified 
by the neutrino chemical potentials /ij^^ (a = e, ii, t) or, by the degeneracy param- 
eters, the ratios of the neutral lepton chemical potentials to the temperature (in 
energy units) = /i^^/fcT, where 



Prior to annihilation, Tiy = T^, while post-e^ annihilation {T,y/Tj)^ = A/11. 
Although in principle the asymmetry among the different neutrino flavors may be 
different, mixing among the three active neutrinos {u,., v^, Vt) ensures that at BBN, 
ie ~ -t-Ai ~ Lt (^e ~ ^/ij ~ ^r)^- If L^, is measured post-e^ annihilation, as is t/b, 
then for < 1, « SL^^ and, for ^ = < 1, « 0.75^. 

Although any neutrino degeneracy (^i^^ < as well as > 0) increases the energy 
density in the relativistic neutrinos, resulting in an effective AN^, ^ (see eq. 10), 
the range of |^| of interest to BBN is limited to sufficiently small values that the 
increase in S due to a non-zero ^ is negligible. However, a small asymmetry between 
electron type neutrinos and antineutrinos (^e ^ 10~^; L ^ 0.007), while large com- 
pared to the baryon asymmetry, can have a significant impact on BBN since the 
affect the interconversion of neutrons to protons. A non-zero results in different 
(compared to SBBN) numbers of Ue and 9^, altering the n/p ratio at BBN, thereby 
changing the yields (compared to SBBN) of the light nuclides. 

Of the light, relic nuclei, the neutron limited ''He abundance is most sensitive 
to a non-zero ^He is a good "leptometer" . In concert with the abundances of D, 
^He, and ""Li, which are good baryometers, the ^He abundance provides a test of the 
consistency of the standard model along with constraints on non-standard models. 
The analytic fits presented below (§2.2) are reasonably accurate for in the range, 
—0.1 ~ 0.1, corresponding to a total lepton number limited to \L\ < 0.07. 

While this may seem small, recall that a similar measure of the baryon asymmetry 
is orders of magnitude smaller: t/b ~ 6 x 10~^°. 




(12) 
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2. An Overview of Primordial Nucleosynthesis 

The early 1 ms), hot, dense Universe is filled with radiation (7s, pairs, i^s of all 
flavors), along with dynamically and numerically insignificant amounts of baryons 
(nucleons) and dark matter particles. Nuclear and weak interactions are occurring 
among the neutrons, protons, , and ps {e.g., n+p < — > D + ^] p + e~ < — > n + Vg) 
at rates fast compared to the universal expansion rate. At such high temperatures 
(T ^ 3 MeV), in an environment where the nucleon to photon ratio is very small 
(?7io 3 — 10), the abundances of complex nuclei (D, ^He,^He, ^Li) are tiny in 
comparison to those of the free miclcons (neutrons and protons). At the same time, 
the charged-current weak interactions are regulating the neutron to proton ratio, 
initially keeping it close to its equilibrium value 

{n/p%, = e-^™/^, (13) 

where Am is the neutron - proton mass (energy) difference. In this context it is 
worth noting that if there is an asymmetry between the numbers of and 
the equilibrium neutron-to-proton ratio is modified to {n/p)eq = exp(— Am/T — 
Me/T) = e-««(n/p)0,. 

As the Universe expands and cools, the lighter protons are favored over the heav- 
ier neutrons and the neutron-to-proton ratio decreases, tracking the equilibrium 
form in eq. 13. But, as the temperature decreases below T ~ 0.8 MeV, when the Uni- 
verse is ~ 1 second old, the weak interactions are too slow to maintain equilibrium 
and the neutron-to- proton ratio, while continuing to fall, deviates from [exceeds) 
the equilibrium value. Since the n/p ratio depends on the competition between the 
weak interaction rates and the early-Universe expansion rate (as well as on a possi- 
ble neutrino asymmetry), deviations from the standard model {e.g., pn.^ Pk + Px 
or 7^ 0) will change the relative numbers of neutrons and protons available for 
building the complex nuclides. 

As noted above, while neutrons and protons are interconverting, they are also 
colliding among themselves creating complex nuclides, e.g., dcutorons. However, at 
early times, when the density and average energy of the CBR photons are very high, 
the newly formed deuterons find themselves bathed in a background of high-energy 
gamma rays capable of photodissociating them. Since there are more than a billion 
CBR photons for every nucleon in the Universe, the deuteron is photodissociated 
before it can capture a neutron (or a proton, or another deuteron) to build the 
heavier nuclides. This bottleneck to BBN persists until the temperature drops suffi- 
ciently below the binding energy of the deuteron, when there are too few photons 
energetic enough to photodissociate them before they capture miclcons, launching 
BBN. This transition (smooth, but rapid) occurs after annihilation, when the 
Universe is a few minutes old and the temperature has dropped below ~ 80 keV. 

Once BBN begins in earnest, neutrons and protons quickly combine to build 
D, ^H, '^He, and ''He. Since there are no stable mass-5 nuclides, a new bottleneck 
appears at ^He. Nuclear reactions quickly incorporate all available neutrons into 
^He, the most strongly bound of the light nuclides. Jumping the gap at mass-5 
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Fig. 1. The SBBN-predicted primordial abundances of D, ^He, ^lA (relative to hydrogen by num- 
ber), and the "^He mass fraction (Yp), as functions of the baryon abundance parameter ??io. The 
widths of the bands reflect the uncertainties in the nuclear and weak interaction rates. 

requires Coulomb suppressed reactions of ''He with D, or ^H, or ^He, guaranteeing 
that the abundances of the heavier nuchdes are severely depressed below that of 
^He (and even of D and '^He), and that the ^He abundance is determined by the 
neutron abundance when BBN begins. The few reactions that manage to bridge the 
mass-5 gap lead mainly to mass-7 (^Li or, to ^Be which, later, when the Universe 
has cooled further, will capture an electron and decay to ^Li); for the range of 77B of 
interest, the BBN-predicted abundance of ^Li is more than 3 orders of magnitude 
below that of the more tightly bound ^Li. Finally, there is another gap at mass-8, 
ensuring that there is no astrophysically significant production of heavier nuclides. 

The primordial nuclear reactor is short-lived. As the temperature drops below 
T ^ 30 keV, when the Universe is ~ 20 minutes old, Coulomb barriers abruptly 
suppress all nuclear reactions. Afterwards, until the first stars form, no pre-existing, 
primordial nuclides are destroyed (except for those like and ^Be that are unstable 
and decay) and no new nuclides are created. In ^ 1000 seconds BBN has run its 
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course. 

With this as background, the trends of the SBBN-predicted primordial abun- 
dances of the Hght nucUdes with baryon abundance shown in Figure 1 can be 
understood. The reactions burning D and ^He (along with ^H) to ^He are very 
fast (compared to the universal expansion rate) once the deuterium bottleneck is 
breached, ensuring that almost all neutrons present at that time are incorporated 
into *He. As a result, since ^He production is not rate limited, its primordial abun- 
dance is very insensitive (only logarithmically) to the baryon abundance. The very 
slight increase in Yp with increasing rjB reflects the fact that for a higher baryon 
abundance BBN begins slightly earlier, when slightly more neutrons are available. 
The thickness of the Yp curve in Fig. 1 reflects the very small uncertainty in the 
BBN prediction; the uncertainty in Yp (^ 0.2%; ay ^ 0.0005) is dominated by the 
very small error in the weak interaction rates which are normalized by the neu- 
tron lifetime (r„ = 885.7 ± 0.8 s). The differences among the Yp predictions from 
independent BBN codes are typically no larger than AYp ^ 0.0002. 

Nuclear reactions burn D, ^H, and '^He to ^He, the most tightly bound of the 
light nuclides, at a rate which increases with increasing nucleon density, accounting 
for the decrease in the abundances of D and '^He (the latter receives a contribution 
from the /3-decay of ^H) with higher values of rjB- The behavior of ''Li is more inter- 
esting, reflecting two pathways to mass-7. At the relatively low values of r]io ^ 3, 
mass-7 is largely synthesized as ''Li by ^H(a,7)''Li reactions. ''Li is easily destroyed 
in collisions with protons. So, for low nucleon abimdance, as rjB increases, destruc- 
tion is faster than production and ''Li/H decreases. In contrast, at relatively high 
values of r]io ^ 3, mass-7 is largely synthesized as ''Be via ^He(a,7)''Be reactions. 
^Be is more tightly bound than ''Li and, therefore, harder to destroy. As tib increases 
at high nucleon abundance, the primordial abundance of ''Be increases. Later in the 
evolution of the Universe, when it is cooler and neutral atoms begin to form, '^Be 
captures an electron and /3-decays to ''Li. These two paths to mass-7 account for 
the valley shape of the ''Li abundance curve in Fig. 1. 

Not shown on Figure 1 are the BBN-predicted relic abundances of ^Li, ^Be, ^°B, 
and ^^B. Their production is suppressed by the gap at mass-8. For the same range 
in r]B, all of them lie offscale, in the range 10~^° — 10~^^. 

For SBBN the relic abundances of the light nuclides depend on only one free 
parameter, the nucleon abundance parameter rjB- As Figure 1 reveals, for the "in- 
teresting" range (see below) of 4 ^ 7710 ^ 8, the ^He mass fraction is expected to be 
Yp w 0.25, with negligible dependence on r/e while D/H and ''He/H decrease from 
w 10"'' to w IQ-s, and ''Li/H increases from « 10"^" to w 10"". The light rniclide 
relic abundances span some nine orders of magnitude, yet if SBBN is correct, one 
choice of jjb (within the errors) should yield predictions consistent with observa- 
tions. Before confronting the theory with data, it is useful to consider a few generic 
examples of BBN in the presence of nonstandard physics and/or cosmology. 
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2.1. Nonstandard BBN 

The variety of modifications to the standard models of particle physics and of cos- 
mology is very broad, limited only by the creativity of theorists. Many nonstandard 
models introduce several, new, free parameters in addition to the baryon abundance 
parameter rjB- Since there arc only four nuclides whose relic abundance is large 
enough to be astrophysically interesting and, as will be explained below in more 
detail, only three for which data directly relating to their primordial abundances 
exist at present (D, ^He, ^Li), nonstandard models with two or more additional 
parameters may well be unconstrained by BBN. Furthermore, as discussed in the 
Introduction (see §1.2 and §1.3), there already exist two additional parameters with 
claims to relevance: the expansion rate parameter S (or, AA^,^; see eqs. 8,11) and 
the lepton asymmetry parameter L (or, ^; see eq. 12). 

The primordial abundance of ^Hc depends sensitively on the pre- and the post- 
ed annihilation early universe expansion rate (the Hubble parameter H) and on 
the magnitude of a — Ve. asymmetry because each will affect the n/p ratio at 
BBN (sec, e.g., Stcigman, Schramm & Gunn 1977 (SSG)^; for recent results see 
Kneller & Steigman 2004 (KS)^). A faster expansion {S > 1; AA^,, > 0) leaves less 
time for neutrons to convert into protons and the higher neutron abundance results 
in increased production of ^Hc. For small changes at fixed ijb, AYp « 0.16(5' — 
1) « 0.013 ANi, (KS). Although the relic abimdances of D and ^He do depend 
on the competition between the nuclear reaction rates and the post-e^ annihilation 
expansion rate (faster expansion less D and ^He destruction ^ more D and ^He), 
they arc much less sensitive to relatively small deviations from S* = 1 (ATV^ = 0)®. 
For mass-7 the effect of a nonstandard expansion rate is different at low and high 
values of tjb- At low baryon abundance {ijio < 3), a faster expansion leaves less 
time for ^Li destruction and the relic abundance of mass-7 increases. In contrast, 
at high baryon abundance (ryio Jc^ 3), S* > 1 leaves less time for ^Be production and 
the relic abundance of mass-7 decreases. As for D and ^He, the quantitative change 
in the ''Li abundance is small for small deviations from SBBN. 

For similar reasons, Yp is sensitive to an asymmetry in the electron neutrinos 
which, through the charged current weak interactions, help to regulate the n/p 
ratio. For > 0, there are more neutrinos than antineutrinos, so that reactions 
such as n + Ue — > p + e~, drive down the n/p ratio. For small asymmetry at fixed 
1]B, KS find AYp ?a — 0.23^e. The primordial abimdances of D, •'He. and ''Li, while 
not entirely insensitive to neutrino degeneracy, are much less affected by a nonzero 
than is *He {e.g., Kang & Steigman 1992^°). 

Each of these nonstandard cases (5 ^ 1, ^ ^ 0) will be considered below. While 
certainly not exhaustive of the nonstandard models proposed in the literature, 
they actually have the potential to provide semi-quantitative, if not quantitative, 
understanding of BBN in a large class of nonstandard models. Note that data 
constraining the primordial abundances of at least two different relic nuclei (one 
of which should be ^He) are required to break the degeneracy between the baryon 
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density and the additional parameter resulting from new physics or cosmology. *He 
is a poor baryometer but a very good chronometer and/or, leptometer; D, '^He, ^Li 
have the potential to be good baryometers. 

2.2. Simple - But Accurate - Fits To The Primordial Abundances 

While BBN involves only a limited number of coupled differential equations, they 
are non- linear and not easily solved analytically. As a result, detailed comparisons of 
the theoretical predictions with the inferred relic abundances of the light nuclei re- 
quires numerical calculations, which may obscure key relations between abundances 
and parameters, as well as the underlying physics. In particular, the connection be- 
tween the cosmological parameter set {tib, S, ^e} the abundance data set {yn, 
Yp, yu}^ may be blurred, especially when attempting to formulate a quantitative 
understanding of how the latter constrains the former. However, it is clear from 
Figure 1 that the relic, light nuclide abundances are smoothly varying, monotonic 
functions of t/b over a limited but substantial range. While the BBN-predicted pri- 
mordial abundances are certainly not linearly related to the baryon density (nor 
to the other parameters 5* and ^e), over the restricted ranges identified above, KS^ 
found linear fits to the predicted abundances (or, to powers of them) which work 
very well indeed. Introducing them here enables and simplifies the comparison of 
theory with data (below) and permits a quick, reasonably accurate, back of the 
envelope, identification of the successes of and challenges to BBN. 

For the adopted range of ?7b, J/d = yuiVB) is well fit by a power law, 

yE'^ ^ 46.5r;ro^-«. (14) 

While the true t/d - Vb relation is not precisely a power law, this fit (for 4 < f]io ^ 8) 
is accurate (compared to a numerical calculation) to better than 1%, three times 
smaller than the ~ 3% BBN uncertainty estimated by Buries, Nollett, Turner 2001 
(BNT)^-'^; this fit and the numerical calculation agree with the BNT result to 2% 
or better over the adopted range in 773. Note that since different BBN codes are 
largely independent and often use somewhat different nuclear reaction data sets, the 
differences among their predicted abundances may provide estimates of the overall 
uncertainties. It is convenient to introduce a "deuterium baryon density parameter" 
rju , the value of r/io corresponding to an observationally determined primordial D 
abundance. 

,46.5(1 ±0.03) 
VB = ( 

Generalizing this to include the two other parameters, KS find 

Vn = Vio-Q{S-l) + ^. (16) 

''Yp is the '*He mass fraction while the other abundances are measured by number compared to 
hydrogen. For numerical convenience, yu = 10^(D/H) and j/li = 10^''(Li/H). 
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Fig. 2. Isoabundance curves for Deuterium (dashed lines) and Helium-4 (solid lines) in the expan- 
sion rate factor (S) — baryon abundance (r;io) plane. The *He curves, from bottom to top, are for 
Yp = 0.23, 0.24, 0.25. The D curves, from left to right, are for yu = 4.0, 3.0, 2.0. 



This fit works quite well for 2 ^ j/d ^ 4, corresponding to 5 ^ ?7d 7. In Figure 2 
the deuterium isoabundance curves are shown in the S — r/io plane, while Figure 3 
shows the same isoabundance contours in the ^e^^yio plane. It is clear from Figures 2 
and 3 that D is a sensitive baryometer since, for these ranges of S and Vb ~ Vio- 

Next, consider ^He. While over a much larger range in r/io, Yp varies nearly 
logarthmically with the baryon density parameter, a linear fit to the Yp versus 7710 
relation is actually remarkably accurate over the restricted range considered here. 

Y^^'^ = 0.2384 + O.0016?7io = 0.2384 + ?7io/625. (17) 

Over the same range in 7710 this fit agrees with the numerical calculation and with 
the BNT" predictions for Yp to within 0.0002 (< 0.1%), or better. Any differences 
between this fit and independent, numerical calculations are smaller (much smaller) 
than current estimates of the errors in the observationally inferred primordial value 
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Fig. 3. As in Figure 2, in the neutrino asymmetry (^e) - baryon abundance (J7io) plane. The *He 
curves, from bottom to top, are for Yp = 0.25, 0.24, 0.23. The D curves, from left to right, are 
for yu = 4.0, 3.0, 2.0. 



of Yp. The following linear fits, including the total error estimate, to the Yp - S 
and Yp - relations from KS work very well over the adopted parameter ranges 
(see Figures 2 & 3). 

Y^^^ = 0.2384 ± 0.0006 + O.OOlGr/io + 0.16(5' - 1) - 0.23^6- (18) 

As an aside, the dependence of the "^He mass fraction on the neutron lifetime (r„) 
can be included in eq. 18 by adding a term 0.0002(r„ — 887.5), where t„ is in 
seconds. A very recent, new measurement of t„ by Serebrov et al}'^ suggests that the 
currently accepted value (r„ = 887.5 s) should be reduced by 7.2 s. If confirmed, this 
would lead to a slightly smaller BBN-predicted ''He abundance: AYp ~ —0.0014. 
The corresponding shift in the ^He inferred baryon density parameter is negligible 
compared to its range of uncertainty (A77B/77B — —0.14), as is that for the shift in 
the upper bound to Ni, (AN™"^ — +0.11). These corrections are ignored here. 
In analogy with the deuterium baryon density parameter introduced above, it 
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is convenient to introduce r^He, defined by 

r/He = 625(Yp - 0.2384 ± 0.0006), (19) 

so that 

rme = mo + 100(5-1)-^. (20) 

For SBBN (5 = 1 & = 0)) ?7He is the value of ?7io corresponding to the adopted 
value of Yp. Once Yp is chosen, the resulting vahic of t/hc provides a linear con- 
straint on the combination of ryio, S, and in eq. 20. This fit works well^ for 
0.23 < Yp < 0.25, corresponding to -5 < r?He $ 7. As Figures 2 & 3 reveal, *He 
is an excellent chronometer and/or leptometer, since the Yp isoabundance curves 
are nearly horizontal (and very nearly orthogonal to the deuterium isoabundance 
curves) . 

As with D, the ^Li abundance'^ is well described by a power law in ryio over the 

range in baryon abundance explored here: j/li ^ 10"'''(Li/H) oc ry^Q. The following 
KS fit agrees with the BBN predictions to better than 3% over the adopted range 
in mo, 

vir ^ %■ (21) 

While this fit predicts slightly smaller lithium abundances compared to those of 
BNT^^, the differences are at the 5-8% level, small compared to the BNT uncer- 
tainty estimates as well as those of Hata et al. (1995)^^ (~ 10 — 20%). 

In analogy with ?7d and ?7ho defined above, the lithium baryon abundance pa- 
rameter 77Li (allowing for a 10% overall uncertainty) is defined by 

r?Li = (8.5(1 ±0.1)yLi)'/'. (22) 

The simple, linear relation for i]u as a function of 7710, S, ^e, which KS find fits 
reasonably well over the adopted parameter ranges is, 

Vu = mo-S{S-l)-^. (23) 

This fit works well for 3 ^ yLi ^ 5, corresponding to 5 ^ 77Li ^ 7, but it breaks down 
for 2/Li ~ 2 {rju ^ 4); see Fig. 1. As is the case for deuterium, lithium can be an 
excellent baryometer since, for the restricted ranges of S and under consideration 
here, rju « ?7io- 

Finally, it may be of interest to note that for ^He the power law j/3 — rjB relation, 
where 2/3 = 10^(^He/H), which is reasonably accurate for 4 < r]io ^ 8 is 

j/3 = 3.1(1 ± 0.03)77ro°-6. (24) 

The difficulty of using current observational data, limited to chemically evolved 
regions of the Galaxy, to infer the primordial abundance of ^He, along with the 



"^It is common in the astronomical literature to present the hthium abundance logarithmically: 
[Li] s 12+Iog(Li/H) = 2 + log(yLi). 
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relatively weak dependence of on 7710, limits the utility of this nuclide as a 
baryometer^"'. "^He can, however, be used as a test of BBN consistency. 

2.3. SBBN-Predicted Primordial Abundances 

Before discussing the current status of the observationally determined abundances 
(and their uncertainties) of the light nuclides, it is interesting to assume SBBN 
and, for the one free parameter, rj-Q, use the value inferred from non-BBN data 
such as the CBR (WMAP) and Large Scale Structure (LSS)^ to predict the relic 
abundances. 

From WMAP alone, Spergel et al. 2003^ derive 7710 = 6.3 ± 0.3. Using the 
fits from §2.2, with S = 1 and = 1) the SBBN-predicted relic abundances are: 
yo = 2.45 ± 0.20; 2/3 = 1-03 ± 0.04; Yp = 0.2485 ± 0.0008; = 4.67 ± 0.64 
([Li]p = 2.67 ±0.06). 

When Spergel et al. 2003^ combine the WMAP CBR data with those from Large 
Scale Structure, they derive a consistent, but slightly smaller (slightly more precise) 
baryon abundance parameter 7710 = 6.14±0.25. For this choice the SBBN-predicted 
relic abundances are: yo = 2.56 ± 0.18; = 1.04 ± 0.04; Yp = 0.2482 ± 0.0007; 
2/Li = 4.44 ± 0.57 ([Li]p = 2.651°;°^). 

3. Observationally Inferred Primordial Abundances 
3.1. The Primordial Deuterium Abundance 

Deuterium is the baryometer of choice since its post-BBN evolution is simple (and 
monotonic!) and its BBN-predicted relic abundance depends sensitively on the 

baryon abundance (r/D oc ?7g^'^). As the most weakly bound of the light nuclides, any 
deuterium cycled through stars is burned to "^He and beyond during the pre-main 
sequence, convective (fully mixed) evolutionary stage^^. Thus, deuterium observed 
anywhere, anytime, should provide a lower bound to the primordial D abundance. 
For "young" systems at high redshift and/or with very low metallicity, which have 
experienced very limited stellar evolution, the observed D abundance should be 
close to the primordial value. Thus, although there are observations of deuterium 
in the solar system and the interstellar medium (ISM) of the Galaxy which pro- 
vide interesting lower bounds to the primordial abundance, it is the observations 
of relic D in a few (too few!), high redshift, low metallicity, QSO absorption line 
systems (QSOALS) which are of most value in enabling estimates of its primordial 
abundance. 

In contrast to the great asset of the simple post-BBN evolution, the identical 
absorption spectra of Di and Hi (modulo the velocity/wavelength shift resulting 
from the heavier reduced mass of the deuterium atom) is a severe liability, limiting 
drastically the number of useful targets in the vast Lyman-alpha forest of QSO 
absorption spectra (see Kirkman et al}^ for further discussion). As a result, it is 
essential to choose target QSOALS whose velocity structure is "simple" since a 
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Fig. 4. The observationally inferred primordial deuterium abundances (the ratio of D to H by 
number) versus a logarithmic measure of the metallicity, relative to solar ([Si/H]), for five high 
redshift, low metallicity QSOALS (filled circles) through 2003. The error bars are the quoted Icr 
uncertainties. Also shown for comparison are the D/H ratios inferred from observations of the local 
interstellar medium (ISM; filled square) and that for the pre-solar nebula (Sun; filled triangle). 



low column density Hi absorber, shifted by ~ 81 km/s with respect to the main 
Hi absorber (an "interloper") could masquerade as Di absorption^^. If this de- 
generacy is not recognized, a D/H ratio which is too high could be inferred. Since 
there are many more low-column density absorbers than those with high H i column 
densities, absorption systems with somewhat lower Hi column density {e.g., Lyman- 
limit systems: LLS) may be more susceptible to this contamination than the higher 
Hi column density absorbers {e.g., damped Lya absorbers: DLA). While the DLA 
do have many advantages over the LLS, a precise determination of the Hi col- 
umn density utilizing the damping wings of the H i absorption requires an accurate 
placement of the continuum, which could be compromised by Hi interlopers. This 
might lead to errors in the Hi column density. These complications are real and 
the path to primordial D using QSOALS has not been straightforward, with some 
abundance claims having had to be withdrawn or revised. Through 2003 there were 
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Fig. 5. Figure 4 updated to 2004 to include the one new deuterium abundance determination^® 
for a high redshift, low metallicity QSOALS. The dashed lines show the SBBN-predicted Icr band 
for the WMAP baryon abundance. 



only five "simple" QSOALS with deuterium detections leading to reasonably robust 
abundance determinations^^ (and references therein); these are shown in Figure 4 
along with the corresponding solar system and ISM D abundances. It is clear from 
Figure 4, that there is significant dispersion among the derived D abundances at 
low metallicity which, so far, masks the anticipated primordial deuterium plateau, 
suggesting that systematic (or random) errors, whose magnitudes are hard to es- 
timate, may have contaminated the determinations of the Di and/or Hi column 
densities. 

It might be hoped that as more data are acquired, the excessive dispersion 
among the deuterium abundances seen in Fig. 4 would decrease leading to a better 
defined deuterium plateau. Be careful what you wish for! In 2004, Crighton et al}^ 
identified deuterium absorption in another high redshift, low metallicity QSOALS 
and derived its abundance. The updated version of Fig. 4 is shown in Fig. 5. Now, 
the dispersion is larger! The low value of the Crighton et al}^ D abundance. 
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similar to that in the more highly evolved ISM and the pre-solar nebula, is puzzling. 
D. Tytler (private communication) and colleagues have data (unpublished) for the 
same QSOALS and they find acceptable fits with lower Hi column densities and 
no visible D i; perhaps this system is contaminated by an inter loper-'^^. The two low 
D/H ratios, if not affected by random or systematic errors, may be artifacts of a 
small statistical sample or, they may have resulted from "young" regions in which 
some of the relic deuterium has been destroyed or depleted onto dust^^. These 
suggestions pose challenges since any cycling of gas through stars should have led 
to an increase in the heavy element abundances (these QSOALS are metal poor) 
and, at low metallicity the amount of dust is expected to be small. 

For the Spergel et al? baryon abundance of r^io = 6.14 ± 0.25, the SBBN- 
predicted deuterium abundance is j/d = 2.6 ± 0.2. As may be seen in Fig. 5, the 
current data exhibit the Goldilocks effect: two D/H ratios are "too small", two are 
"too large", and two are "just right". Nonetheless, it is clear that the sparse data 
currently available are in very good agreement with the SBBN prediction. If the 
weighted mean of D/H for the six QSOALS is adopted, but the dispersion in the 
mean is used in place of the error in the mean, = 2.4 ± 0.4, corresponding to 
an SBBN-predictcd baryon abundance of 7710 — 6.4 ± 0.6, in excellent agreement 
with the WMAP value. Were it not for the excessive dispersion among the extant 
deuterium abundance determinations, the precision of the baryon abundance de- 
termined from SBBN would be considerably enhanced. For example, if the formal 
error in the mean is used, 7710 = 6.4 ± 0.2 is the deuterium-based, SBBN predic- 
tion, very close to that for the baryon abundance derived from the GBR alone^ 
{r]\o^^^ = 6.3 ± 0.3). More data is crucial to deuterium fulfilling its potential. 

3.2. The Primordial Helium- 3 Abundance 

The post-BBN evolution of ^He, involving competition among stellar production, 
destruction, and survival, is considerably more complex and model dependent than 
that of D. Interstellar ^He incorporated into stars is burned to ''He (and beyond) in 
the hotter interiors, but it is preserved in the cooler, outer layers. Furthermore, while 
hydrogen burning in cooler, low-mass stars is a net producer of ^He^°, it is unclear 
how much of this newly synthesized '^Hc is returned to the interstellar medium and 
how much of it is consumed in post-main sequence evolution {e.g., Sackmann & 
Boothroyd^^). For years it had been anticipated that net stellar production would 
prevail in this competition, so that the ^He abundance would increase with time 
(and with metallicity)^^. 

Observations of '^He, are restricted to the solar system and the Galaxy. Since 
for the latter there is a clear gradient of metallicity with location, a gradient in 
"^Hc abundance was also expected. However, as is clear from Figures 6 & 7, the 
data^^'^^ reveal no statistically significant correlation between the ^He abundance 
and metallicity or location in the Galaxy, suggesting a very delicate balance between 
net production and net destruction of ^He. For a recent review of the current status 
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Fig. 6. The '^He abundance determinations (by number relative to H) in the ISM of the Galaxy 
(from BRB^*) as a function of the corresponding oxygen abundances. The solar symbol indicates 
the ^He abundance for the pre-solar nebula. The dashed lines show the Itr band adopted by BRB. 



of '^He evolution, see Romano et al.^^. 

While the absence of a gradient suggests the mean ("plateau") ^He abundance 
in the Galaxy (7/3 « 1.9 ± 0.6) might provide a good estimate of the primordial 
abundance, Bania, Rood & Balser (BRB)^'* prefer to adopt as an upper limit to 
the primordial abundance, the ^He abundance measured in the most distant (from 
the Galactic center), most metal poor, Galactic Hii region, 7/3 1.1 ±0.2; see Figs. 
6 & 7. This choice is in excellent agreement with the SBBN/WMAP predicted 
abundance of ya = 1.04 ± 0.04 (see §2.2). While both D and ^He are consistent 
with the SBBN predictions, '^He is a less sensitive baryometer than is D since 
(D/H)bbn oc 773^-^, while (3He/H)BBN oc 773" ^ For example, if = 1.1 ± 0.2 is 
adopted for the "^He primordial abundance, 77io('^He) = 6.0 ± 1.7. While the central 
value of the ^He-inferred baryon density parameter is in nearly perfect agreement 
with the WMAP value^ , the allowed range of rfs is far too large to be very useful. 
Still, •^He can provide a valuable BBN consistency check. 
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Fig. 7. As in Figure 6 but, for the '^He abundances as a function of distance from the center of 
the Galaxy. 



3.3. The Primordial Helium-^ Abundance 

The post-BBN evolution of *He is quite simple. As gas cycles through generations 
of stars, hydrogen is burned to helium-4 (and beyond), increasing the "^He abun- 
dance above its primordial value. The "^He mass fraction in the Universe at the 
present epoch, Yq, has received a significant contribution from post-BBN, stellar 
nucleosynthesis, so that Yq > Yp. However, since the "metals" such as oxygen are 
produced by short-lived, massive stars and "^He is synthesized (to a greater or lesser 
extent) by all stars, at very low metallicity the increase in Y should lag that in 
e.g., 0/H so that as 0/H ^ 0, Y ^ Yp. As is the case for deuterium and lithium, 
a '^He "plateau" is expected at sufficiently low metallicity. Therefore, although ''He 
is observed in the Sun and in Galactic H ii regions, the key data for inferring its pri- 
mordial abundance are provided by observations of helium and hydrogen emission 
(recombination) lines from low-metallicity, extragalactic Hii regions. The present 
inventory of such regions studied for their helium content exceeds 80 (see Izotov & 
Thuan (IT)^^). Since with such a large data set even modest observational errors 
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Fig. 8. The observationally inferred primordial *He mass fractions from 1978 until 2004. The error 
bars are the quoted Itr uncertainties. Also shown is the SBBN-predicted relic abundance (solid 
line) for the WMAP baryon abundance, along with the la uncertainty (dashed lines) of the SBBN 
prediction. 



for the individual H ii regions can lead to an inferred primordial abundance whose 
formal statistical uncertainty may be quite small, special care must be taken to 
include hitherto ignored or unaccounted for systematic corrections and/or errors. 
It is the general consensus that the present uncertainty in Yp is dominated by 
the latter, rather than by the former errors. Indeed, many of the potential pitfalls 
were identified by Davidson & Kinman^^ in a prescient, 1985 paper. In the abstract 
they say, "The most often quoted estimates of the primordial helium abundance 
are optimistic in the sense that quoted uncertainties usually do not include some 
potentially serious systematic errors." 

To provide a context for the discussion of the most recent data and analyses. 
Figure 8 offers a compilation of the history of Yp determinations^^'^^ derived us- 
ing data from low metallicity, extragalactic Hii regions. Notice that all of these 
estimates, taken at face value, fall below the SBBN/WMAP predicted primordial 
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abundance by at least 2(t, recmphasizing the importance of accounting for system- 
atic uncertainties. With this in mind, we turn to recent reanalyses^^'"^" of the IT 
data^^, supplemented by key observations of a local, higher metallicity H ii region^°. 

Prior and subsequent to the Davidson & Kinman papcr^^ astronomers have 
generally been aware of the important sources of potential systematic errors asso- 
ciated with using recombination line data to infer the helium abundance. However, 
attempts to account for them have often been unsystematic or, entirely absent. The 
current conventional wisdom that the accuracy of the data demands their inclusion 
has led to some attempts to account for a few of them or, for combinations of a 
few of them26.29,30,3i,32 rphe Olive & Skillman (OS)^^ analysis of the IT data is 
the most systematic to date. Following criteria outlined in their 2001 paper^®, OS 
found they were able to apply their analysis to only 7 of the 82 IT Hii regions. 
This tiny data set, combined with its limited range in metallicity (oxygen abun- 
dance), severely limits the statistical significance of any conclusions OS can extract 
from it. In Figure 9 arc shown the differences between the OS-inferred and the IT- 
inferred helium abundances. For these seven Hii regions there is no evidence that 
AY = Y*-*^ — Y'"^ is correlated with metallicity. The weighted mean offset along 
with the error in the mean arc AY = 0.0029 ± 0.0032 (the average offset and the 
average error are AY = 0.0009 ± 0.0095), consistent with zero at la. 

If the weighted mean offset is applied to the IT-derived primordial abundance 
of Yp"^ = 0.2443 ± 0.0015, the "corrected" primordial value becomes 

Y^^ = Y^^ + AY = 0.2472 ± 0.0035, (25) 

leading to a 2a upper bound on the primordial abundance of Yp^ < 0.254. In 
contrast, OS prefer to fit these seven data points to a linear Y versus 0/H relation 
and, from it, derive the primordial abundance. Their IT-revised abundances, along 
with, for comparison, that from their reanalysis of the Peimbert et al.^^ data for 
an Hii region in the SMC (to be discussed next), are shown in Figure 10. It is not 
surprising that for only seven data points, each with larger errors than those adopted 
by IT, spanning such a narrow range in metallicity, their linear fit. Yy*^ = 0.2495 ± 
0.0092-1- (54± 187)(0/H), is not statistically significant. Indeed, it is not preferred 
over the simple weighted mean of the seven helium abundances (0.252 ± 0.003), 
since the per degree of freedom is actually higher for the linear fit. In fact, there 
is no statistically significant correlation between Y and 0/H for the IT-derived 
abundances for these seven Hii regions either. As valuable as is their reanalysis of 
the IT data, the OS conclusion that Yp = 0.249 ± 0.009 is not supported by the 
sparse data set they used^. Unless and until an analysis is performed of a much 
larger data set, with a longer metallicity baseline, the estimate in eq. 25, and its 
corresponding 2a upper bound, may provide a good starting point at present for 
an approach to the primordial abundance of ^He. 

"^Note that OS used the corresponding la upper bound of 0.258 for the upper bound to their 
"favored" primordial abundance range. 
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Fig. 9. The differences between the OS and IT *He abundances, AY=Y°S_yIT for the OS- 
selected IT H 11 regions versus the corresponding oxygen abundances. The solid line is the weighted 
mean of the helium mass fraction differences, while the dashed line shows the unweighted average 
of the differences. 



Another correction, not directly constrained by the analysis of OS, is related 
to the inhomogeneous nature of Hii regions. Unlike classical, textbook, homoge- 
neous, Stromgren spheres, real Hii regions are filamentary and inhomogeneous, 
with variations in electron density and temperature likely produced by shocks and 
winds from pockets of hot, young stars. Esteban and Peimbert^'^ have noted that 
temperature fluctuations can have a direct effect on the helium abundance derived 
from recombination lines. This effect was investigated theoretically using models of 
Hii regions^^, but more directly by Peimbert et al.^^ using data from a nearby, spa- 
tially resolved, Hii region in the Small Magellanic Cloud (SMC), along with their 
reanalyses of four Hii regions selected from IT. While the SMC Hii region formed 
out of chemically evolved gas and, therefore, cannot be used by itself to derive pri- 
mordial abundances, the spatial resolution it offers permits a direct investigation of 
many potential systematic effects. In particular, since recombination lines are used, 
the observations are blind to any neutral helium or hydrogen. Estimates of the "ion- 
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Fig. 10. The OS-revised ''He versus oxygen abundances for the seven IT H II regions and the SMC 
H II region from PPR. The solid line is the weighted mean of the helium abundances for all eight 
of the Hll regions reanalyzed by OS. 

ization correction factor" {icf), while model dependent, are large"^^. For example, 
using models of H ii regions ionized by distributions of stars of different masses and 
ages and comparing to the IT (1998) data, Gruenwald et al.^^ concluded that IT 
overestimated the primordial ''He abundance by AY'^^^(zc/) « 0.006±0.002; Sauer 
& Jedamzik'^^ find a similar, even larger, correction. If this correction is applied to 
the OS-revised, IT primordial abundance in eq. 25, the new, ic/-corrected value is 

Y^s^(ic/) = Y^^ - AY°s^(ic/) = 0.241 ± 0.004. (26) 

In addition to the subset of 7 of the 82 IT H ii regions which meet their criteria, 
OS also reanalyzed the Peimbert et al. data'^'^ for the SMC Hii region. This OS- 
revised data point is shown in Figure 10 at the highest oxygen abundance. Notice 
that the eight data points plotted in Fig. 10 show no evidence of the expected 
increase of Y with metallicity''; this is likely due to the small sample size. The 

"^Absence of evidence is NOT evidence of absence. 
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weighted mean He abundance for these eight H ii regions is = 0.250 ± 0.002, 
corresponding to a 2a upper bound of Yp < 0.254. Coincidentally, this is the same 
2a upper bound as that found from the mean of the seven IT Hii regions (see 
eq. 25) and, also, the 2a upper bound from the OS-reanalyzed SMC H ii region 
alone. If the ionization corrections from Gruenwald et al.'^^ are applied to each of 
these eight Hii regions, it is found that the mean AY{icf)s = —0.002 ± 0.002, so 
that including this correction, while accounting for the increased error, leaves the 
2(7 upper boimd of Yp < 0.254 unchanged. 

The lesson from the discussion above is that while recent attempts to determine 
the primordial abundance of ^He may have achieved high precision, their accuracy 
remains in question. The latter is limited by our understanding of and our ability 
to account for systematic errors and biases, not by the statistical uncertainties. The 
good news is that carefully organized, detailed studies of only a few (~ a dozen?) 
low metallicity, extragalactic Hii regions may go a long way towards an accurate 
determination of Yp. The bad news is that many astronomers and telescope allo- 
cation committee members are unaware that this is an interesting and important 
problem, worth their effort and telescope time. At present then, the best that can 
be done is to adopt a defensible value for Yp and, especially, its uncertainty. To this 
end, in the following the estimate in eq. 26 is chosen: Yp = 0.241 ±0.004. While the 
central value of Yp is low, it is within 2a (~ 1.75c7) of the SBBN/WMAP expected 
central value of Yp = 0.248 (see §2.2). Note that the extrapolation of the linear 
fit of the {Y, 0/H} data from the lowest mctalhcity (O/H 2 x 10~^) to zero 
metallicity (Yp) corresponds to AY « 0.0009, well within the uncertainties of Yp. 

In setting contraints on new physics, an upper bound to Yp is required. A robust 
upper bound suggested by the above discussion is Yp < 0.254. As an example, the 
SBBN/WMAP lower bound (at - 2a) to Yp is 0.247, so that AYp < 0.007. This 
corresponds to the robust upper bounds S < 1.04 and N,^ < 3.5, eliminating (just 
barely) even one, new, light scalar, and bounding the lepton asymmetry from below: 
> -0.03. 

3.4. The Primordial Lithium-7 Abundance 

In the post-BBN universe ^Li, along with ^Li, ''Be, ^°B, and ^^B, is produced in 
the Galaxy by cosmic ray spallation/fusion reactions. Furthermore, observations 
of super-lithium rich red giants provide evidence that (at least some) stars are 
net producers of lithium. Therefore, even though lithium is easily destroyed in the 
hot interiors of stars, theoretical expectations supported by the observational data 
shown in Figure 11 suggest that while lithium may have been depleted in many 
stars, the overall trend is that its abundance has increased with time. Therefore, 
in order to probe the BBN yield of ^Li, it is necessary to restrict attention to the 
oldest, most metal-poor halo stars in the Galaxy (the "Spite Plateau") seen at 
low metallicity in Fig. 11. Using a selected set of the lowest metallicity halo stars, 
Ryan et al.^^ claim evidence for a 0.3 dex increase in the lithium abundance ([Li] = 
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Fig. 11. Lithium abundances, log e(Li) s [Li] = 12-|-log(Li/H) versus metallicity (on a log scale 
relative to solar) from a compilation of stellar observations by V. V. Smith. The solid line is 
intended to guide the eye to the "Spite Plateau" . 



12+log(Li/H)) for —3.5 < [Fe/H] < —1, and they derive a primordial abundance of 
[Li]p « 2.0—2.1. This abundance is low compared to the value found by Thorburn^^, 
who derived [Li]p « 2.25 ± 0.10. The stellar temperature scale plays a key role in 
using the observed equivalent widths to derive the ''Li abundance. Studies of halo 
and Galactic Globular Cluster stars employing the infrared flux method eflfective 
temperature scale suggest a higher lithium plateau abundance^^: [Li]p = 2.24±0.01, 
similar to Thorburn's'^^ value. Recently, Melendez & Ramirez'^^ reanalyzed 62 halo 
dwarfs using an improved infrared flux method effective temperature scale. While 
they failed to confirm the [Li] - [Fe/H] correlation claimed by Ryan et al.^'^, they 
suggest an even higher relic lithium abundance: [Li]p = 2.37 ±0.05. A very detailed 
and careful reanalysis of extant observations with great attention to systematic 
uncertainties and the error budget has been done by Charbonnel and Primas^*, who 
find no convincing evidence for a Li trend with metallicity, deriving [Li]p = 2.21 ± 
0.09 for their fidl sample and [Li]p = 2.18 ± 0.07 when they restrict their sample 
to unevolved (dwarf) stars. They suggest the Melendez & Ramirez value should 
be corrected downwards by 0.08 dex to account for different stellar atmosphere 
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models, bringing it into closer agreement with their results. To err on the side of 
conservatism, the lithium abundance of Melendez & Ramirez^^, [Li]p = 2.37 ±0.05, 
which is closer to the SBBN expectation, will be adopted in further comparisons. 

There is tension between the SBBN predicted relic abundance of ^Li ([Li]p = 
2.65to'o6> §2-2) and that derived from recent observational data ([Li]p = 
2.37 ± 0.05). Systematic errors may play a large role confirming or resolving this 
factor of two discrepancy. The role of the stellar temperature scale has already been 
mentioned. Another concern is associated with the temperature structures of the 
atmospheres of these very cool, metal-poor stars. This can be important because 
a large ionization correction is needed since the observed neutral lithium is a mi- 
nor component of the total lithium. Furthermore, since the low metallicity, dwarf, 
halo stars used to constrain primordial lithium are among the oldest in the Galaxy, 
they have had the most time to alter (by dilution and/or destruction) their surface 
lithium abundances, as is seen to be important for many of the higher metallic- 
ity stars shown in Fig. 11. While mixing stellar surface material to the interior 
would destroy or dilute any prestellar lithium, the very small observed dispersion 
among the lithium abundances in the low metallicity halo stars (in contrast to the 
very large spread for the higher metallicity stars) suggests this correction may not 
be large enough (^ 0.1 — 0.2 dex at most) to bridge the gap between theory and 
observation; see, e.g., Pinsonneault et al.^^ and further references therein. 

4. Discussion 

The cosmic nuclear reactor was active for a brief epoch in the early evolution of the 
universe. As the Universe expanded and cooled the nuclear reactor shut down after 

^ 20 minutes, having synthesized in astrophysically interesting abundances only 
the lightest nuclides D, ^He, '^He, and ^Li. For the standard models of cosmology 
and particle physics (SBBN) the relic abundances of these nuclides depend on 
only one adjustable parameter, the baryon abundance parameter rjB (the post- 
ed annihilation value of the baryon (nucleon) to photon ratio). If the standard 
models are the correct description of the physics controlling the evolution of the 
universe, the abundances of the four nuclides should be consistent with a single value 
of rjB and this baryon density parameter should also be consistent with the values 
inferred from the later evolution of the universe {e.g., at present as well as ~ 400 
kyr after BBN, when the relic photons left their imprint on the CBR observed by 
WMAP and other detectors). There are, however, two other particle physics related 
cosmological parameters, the lepton asymmetry parameter and the expansion 
rate parameter S, which can affect the BBN-predicted relic abundances. For SBBN 
it is assumed that = and 5 = 1. Deviations of either or both of these parameters 
from their standard model expected values could signal new physics beyond the 
standard model(s). 

The simplest strategy is to test first the predictions of SBBN. Agreement be- 
tween theory and observations would provide support for the standard models. 
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Disagreements are more difficult to interpret in that while they may be opening a 
window on new physics, they may well be due to unaccounted for systematic errors 
along the path from observations of post-BBN material to the inferred primordial 
abundances. Subject to this latter caveat, the confrontation between theory and 
data can provide useful limits to (some of) the parameters associated with new 
physics which complement those from high precision, terrestrial experiments. In 
the comparisons presented below, the abundances (and their inferred uncertain- 
ties) presented in §3 are adopted and compared to the BBN predictions described 
by the simple fits from §2.2. For ^He the SSBN range in ?7b favored by the adopted 
primordial abundance lies outside the range of validity of the simple fit; for ^He and 
SBBN, the best fit and uncertainty in rj-Q is derived from the more detailed BBN 
calculations. While not all models of new physics proposed in the literature can be 
tested in this manner, this approach does offer the possibility of constraining a large 
subset of them and of providing a useful framework for understanding qualitatively 
how many of the others might affect the BBN predictions. 



4.1. SBBN 

The discussion in §3 identifies a set of primordial abundances. Since these choices 
are certainly subjective and likely to change as more data are acquired, along with 
a better understanding of and accounting for systematic errors, the analytic fits 
presented in §2.2 can be; very useful in relating new conclusions and constraints 
to those presented here. The abundances and nominal Icr uncertainties adopted 
here are: yo = 2.6 ± 0.4, j/3 = 1.1 ± 0.2, Yp = 0.241 ± 0.004, and j/li = 2.34to;25 
([Li]p = 2.37 ± 0.05). Yp < 0.254 is adopted for an upper bound (at ~ 2cr) to 
the primordial ''He mass fraction. The corresponding SBBN values of the baryon 
density parameter are shown in Figure 12, along with that inferred from the CBR 
and observations of Large Scale Structure^ (labelled WMAP). 

As Figure 12 reveals, the adopted relic abundances of D and ^He are consistent 
with the SBBN predictions (770 = 6.1 ±0.6, rjs^^ = 6.0 ± 1.7) and both are in excel- 
lent agreement with the non-BBN value^ (^^wmap = 6.14±0.25). If the most recent 
deuterium abundance determination in a high redshift, low metallicty QSOALS^^ 
is included in estimating the relic D abundance, the mean shifts to a slightly lower 
value [yu = 2.4 ± 0.4), corresponding to a slightly higher estimate for the baryon 
density parameter (r^D = 6.4 ± 0.7), which is still consistent with ^He and with 
WMAP. Were it not for the very large dispersion among the D abundance determi- 
nations (sec §3.1), the formal error in the mean (^ 5%) could have been adopted for 
the uncertainty in j/d, leading to a ~ 3% determination oirju, competitive with that 
from WMAP. Due to the very large observational and evolutionary uncertainties 
associated with ^He, its abundance mainly provides a consistency check at present. 
Since the variations of its predicted relic abundance with S and are similar to 
those for D, ^He will not add new information to that from D in the comparisons 
to be discussed below. 
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Fig. 12. The SBBN values for the early universe (~ 20 minutes) baryon abundance parameter 
•qio inferred from the adopted primordial abundances of D, ''He, ''He, and ''^Li (see §3.1-3.4). Also 
shown is the WMAP-derived CBR and LSS value (~ 400 kyr). 

In addition to the successes of D and ^He, Figure 12 exposes a tension between 
WMAP (and D and ^He) and the adopted primordial abundances of ^He and ''Li. 
The la range determined from ^He is low: 2.2 < ryne < 4-3; however, the 2cr range 
is much larger: 1.7 < r/He < 6.4, encompassing the WMAP-inferred baryon density. 
The ^Li inferred baryon density is also low {■>%{ = 4.5 ± 0.3) and here the adopted 
errors appear to be far too small to bridge the gap to D and WMAP. These tensions 
may be a sign of systematic errors introduced when the observational data is used 
to derive the inferred primordial abundances or, it could be a signal of new physics 
beyond the standard models of cosmology and particle physics. 

4.2. Lithium 

As identified above, the SBBN abundances of D and ^He are in agreement with each 
other and with the non-BBN estimate of the baryon density parameter from Large 
Scale Structure and the CBR. However, while the inferred primordial abundance 
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of ''He is less than 2(T away from the SBBN-predictcd value, that of lithium differs 
from expectations by a factor of ~ 2 (or more). It is unlikely that this conflict 
can be resolved through a non-standard expansion rate {S ^ 1) or a non-zero 
Icpton number (^e 7^ 0). The reason is that in the S — ijb and — Vb planes the 
isoabundance curves for D and ''Li are very nearly parallel (see cqs. 16 & 23 in §2.2 
and Figs. 1 & 2 from Kneller & Steigman^), so that once yo is constrained, there 
is very little freedom to modify i/li- This may be seen by combining eqs. 16 & 23 
to relate 1%! to ijd, 

riLi = VB+mS-i)-U (27) 

Thus, for ?7D ^ 6 and IS" - 1| ^ 0.1, |^e| ^ 0.1, r/Li « ??d ^ 6, so that yu ^ 4 
([Li]p > 2.6). 

Nonetheless, a non-standard physics explanation of the lithium conflict is not 
ruled out. Indeed, there arc models where latc-dccaying, massive particles reinitiate 
BBN, modifying the abundances of the light nuclides produced during the first 20 
minutes. For an extensive, yet likely incomplete list (with apologies) of references, 
see Ref. [36] and further references therein. In such models it is quite possible to 
reduce the original BBN abundance of ''Li to bring it into agreement with the 
value inferred from the observational data^'*'^^'^^'^''. However, it is found that when 
the many new parameters available to these models are adjusted to achieve this 
agreement, the modified relic abundance of ^He is much too large (see, e.g., Ellis, 
Olive, and Vangioni''°). 

The difficulty in reconciling the observed and predicted relic abundances of '^Li 
suggests that the problem may be in the stars. It is not at all unexpected that 
the very old halo stars where lithium is observed will have modified their original 
surface abundances, ''Li in particular (see Pinsonneault et al.^^ and Charbonnel and 
Primas^^ for discussions and many additional references). While there is no dearth 
of physical mechanisms capable of destroying or diluting surface lithium, many of 
which are supported by independent observational data, the challenge has been to 
account for the required depletion (factor of 2 - 3) while maintaining a negligible 
dispersion {'^ 0.1 dex) among the "Spite plateau" lithium abundances. 

Another possibility for reconciling the observed and predicted relic abundances 
of ''Li lies in the nuclear physics. After all, given the estimates of uncertainties in the 
cross sections of the key nuclear reactions leading to the production and destruction 
of mass-7, the BBN-predictcd abundance of ''Li is the most uncertain (~ 10 — 20%) 
of all the light nuclides. Perhaps the conflict between theory and observation is the 
result of an error in the nuclear physics. This possibility was investigated by Cyburt, 
Fields, and Olive^^ who noted that some of the same nuclear reactions of importance 
in BBN, play a role in the standard solar model and are constrained by its success 
in accounting for the observed flux of solar neutrinos. While the uncertainty of a key 
nuclear reaction (•^He(a, 7)'^Be) is large (~ 30%), it is far smaller than the factor of 
~ 3 needed to reconcile the predicted and observationally inferred abundances^^ . 
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Fig. 13. The D and ''He isoabundance curves in the S — »yio plane, as in Fig. 2. The best fit point 
and the error bars correspond to the adopted primordial abundances of D and *He. 

Considering the current state of affairs (no successful resolution based on new 
physics; possible reconciliation based on stellar astrophysics), ^Li is not used below 
where the adopted relic abundances of D and ''He are employed to set constraints 
on S and/or ^e- 

4.3. N on- Standard Expansion Rate: S ^ 1 (^e — 0) 

If the lepton asymmetry is very small, of order the baryon asymmetry, then BBN 
depends on only two free parameters, 77B and S (or N^). Since the primordial 
abundance of D largely probes r/B while that of ''He is most sensitive to S (see 
Fig. 2 and eqs. 16 & 20), for each pair of ?/d and Yp values (within reason) there 
will be a corresponding pair of rfs and S values. For the D and ''He abundances 
adopted above (yo = 2.6±0.4, Yp = 0.241 ±0.004) the best fits for rjB and S, shown 
in Figure 13, are for 7710 = 5.9 ± 0.6 and S = 0.96 ± 0.02; the latter corresponds to 
Njy = 2.5 ± 0.3. These values are completely consistent with those inferred from the 
joint constraints on S and rye from WMAP*^. 
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Fig. 14. The D and ^He isoabundance curves in the — »?io plane, as in Fig. 3. The best fit point 
and the error bars correspond to the adopted primordial abundances of D and *He. 



As expected from the discussion in §2.2, the hthium abundance is largely driven 
by the adopted deuterium abundance and is little affected by the small departure 
from the standard expansion rate. For the above best fit values, ?/Li = 4.3 ± 0.9 
([Li]p = 2.63t°:"^). This class of non-standard models [S ^ 1), while reconciling 
''He with D and with the CBR, is incapable of resolving the lithium conflict. 

4.4. Non-Zero Lepton Number: ^ (S — 1) 

While most popular extensions of the standard model which attempt to account 
for neutrino masses and mixings suggest a universal lepton asymmetry comparable 
in magnitude to the baryon asymmetry (^e ~ ^(t^b) ^ 10^^)^, there is no direct 
evidence that nature has made this choice. Although the CBR is blind to a relatively 

^By charge neutrality the charged lepton excess is equal to the proton excess which constitutes 
^ 87% of the baryon excess. Therefore, any significant lepton asymmetry (^e >> »)b) must be 
hidden in the unobserved relic neutrinos. 
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small lepton asymmetry, BBN provides an indirect probe of it^^. As discussed in 
§1.3 & §2. a lepton asymmetry can change the neutron to proton ratio at BBN, 
modifying the light element yields, especially that of ^He. Assuming 5=1 and 
allowing 7^ 1, BBN now depends on the two adjustable parameters 77B and 
which may be constrained by the primordial abundances of D and ^He. Given the 
strong dependence of Yp on and of yu on rjB, these nuclides offer the most 
leverage. In Figure 14 are shown the D and ^He isoabundance curves {i.e., the fits 
from §2.2) in the S^e—VB plane, along with the best fit point (and its la uncertainties) 
determined by the adopted primordial abundances. The best fit baryon abundance, 
r]w = 6.1 ± 0.6 is virtually identical to the SBBN (and WMAP) value. While the 
best fit lepton asymmetry, = 0.031 ± 0.018, is non-zero, it differs from zero by 
less than 2a (as it should since the adopted value of Yp differs from the SBBN 
expected value by less than 2a). 

As expected from the discussion above for 5 7^ 1 and in §2.2, here, too, the 
lithium abundance is largely driven by the adopted deuterium abundance and 
is little affected by the small lepton asymmetry allowed by D and ^He. For the 
above best fit values, the predicted lithium abundance is virtually identical to the 
SBBN/WMAP and 5 7^ 1 values: yu = 4.3 ± 0.9 ([Li]p = 2.64^^ ?^). A lepton 
asymmetry which reconciles *He with D cannot resolve the lithium conflict. 

4.5. An Example: Alternate Relic Abundances for D and ^He 

It is highly likely that at least some of the tension between D and ''He is due to er- 
rors associated with inferring their primordial abundances from the current observa- 
tional data. As a result, in the future the abundances adopted here may be replaced 
by revised estimates. This is where the simple, analytic fits derived by KS^ and pre- 
sented in §2.2 can be of value to those who lack an in-house BBN code. Provided 
that the revised abundances lie in the ranges 2 i/d ^ 4 and 0.23 ^ Yp ^ 0.25, 
these fits will provide quite accurate, back of the envelope estimates of ?7b, S, ^e, 
and of yu- As an illustration, let's revisit the discussion in §4.3 & §4.4, now adopt- 
ing for D the weighted mean deuterium abundance which results when the most 
recent determination^^ is included, yr, = 2.4 ±0.4 (see §3.1), along with, for ^He, 
the helium abundance derived by applying the OS mean offset to the IT-inferred 
primordial value (see eq. 25) without the «c/-correction, Yp = 0.2472 ±0.0035 (see 
§3.3). These alternate abundances correspond to ijb « 6.5±0.7 and ?7He ~ 5.5±2.2. 

For = 0, the new values for the expansion rate factor and the baryon density 
parameter are 5* = 0.991 ± 0.022 (N,, = 2.9 ± 0.3) and 7?io = 6.4 ± 0.6. While the 
N,^ estimate is entirely consistent with N^ = 3, the corresponding ~ 2a upper bound 
(Nj, ^ 3.5) still excludes even one extra light scalar. The baryon density parameter 
is slightly higher than, but entirely consistent with that inferred from the GBR. As 
anticipated from the previous discussion, the predicted lithium abundance hardly 
changes at all, but it does increase slightly to further exacerbate the conflict with 
the observationally inferred value, yu i=s 4.9 ± 1.1 ([Li]p = 2.69^0 ii)- 
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For S' = 0, the new values for the lepton asymmetry parameter and the baryon 
density parameter are = 0.007±0.016 and r?io = 6.5±0.7. The former is consistent 
with no lepton asymmetry (i.e., with ~ 0{r}^)) and the latter is slightly higher 
than, but still entirely consistent with the baryon density parameter inferred from 
the CBR. As expected, here, too, the predicted lithium abundance increases slightly 
from the already too large SBBN value, ^ 4.9 ± 1.1 ([Li]p = 2.69tHi)- 

4.6. Other Non-Standard Models 

Although the parameterization of BEN in terms of S and explored in the previous 
sections encompasses a large set of non-standard models of cosmology and particle 
physics, it by no means describes all interesting extensions of the standard model. 
As already mentioned, there is a class of models where BBN proceeds normally 
but a second epoch of early universe nucleosynthesis is initiated by the late decay 
of a massive particlc^'^. Despite the fact that such models have many more free 
parameters, such as the mass, abundance, and lifetime of the decaying particle, 
the constraints imposed by the observationally inferred relic abundances of D, ■'He, 
"'He, and ^Li are sufficiently strong to challenge them (see, e.g., Ellis, Olive & 
Vangioni^°). 

There are other models which cannot be simply described by the {tjb, S, ^e} 
parameter set. In most cases they introduce several free parameters in addition 
to the baryon density parameter. Since there are only four nuclides whose relic 
abundances are reasonably constrained, the leverage of BBN on these models may 
be limited in some cases. A case in point is the class of models where the universe is 
inhomogeneous at BBN (IBBN); see the recent article by Lara^^ and the extensive 
references to earlier work therein. In IBBN the geometry of the inhomogeneities 
(spheres, cylinders, ...) is important, as are the scales of the imhomogeneities and 
their amplitudes (density contrasts). Nonetheless, even with all these adjustable 
parameters, except when they take on values indistinguishable from SBBN, IBBN 
models predict an excess of lithium (even more of an excess than for SBBN). This 
is inevitable since in IBBN ^Li is overproduced in the low nucleon density regions 
and ''Be is overproduced in the high density regions (see the multi-valued lithium 
abundance curve in Fig. 1). 

5. Summary and Conclusions 

As the early Universe evolved from hot and dense to cool and dilute, it passed 
through a short-lived epoch when conditions of temperature and density permitted 
the synthesis of astrophysically interesting abundances of D, ^He, ^He, and ''Li. 
At present, observations of these nuclides in a variety of astrophysical sites (stars, 
Galactic and extragalactic Hii regions, QSOALS, etc.) have permitted quite pre- 
cise estimates of their primordial abundances, opening a window on early- universe 
cosmology and providing constraints on physics beyond the standard models of 
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cosmology and particle physics. The relic abundances of D, ^He, and ''Li are nu- 
clear reaction rate limited, providing good probes of the nucleon density when the 
universe was less than a half hour old. This is complementary to baryon density 
estimates from the CBR, whose information was encoded some 400 kyr later, and 
to those provided by observations of the current or recent universe, 13 Gyr after 
BBN. Of these light nuclides, D is the baryometer of choice and current estimates 
of the relic abundance, = 2.6 ±0.4, suggest a baryon to photon ratio, unchanged 
from BBN to the present, of 7710 = 6.1 ±0.6, in excellent agreement with the WMAP 
and LSS determined value of 7710 = 6.14 ± 0.25. For this choice {S = 1, = 0)) 
the less well constrained relic abundance of ^He (2/3 = 1.1 ± 0.2)^^ is also in agree- 
ment with its SBBN expected value. These successes for the standard models are 
tempered by the challenges posed by ^He and ''Li. The SBBN predicted ^He mass 
fraction, Yp = 0.248, differs from the observationally inferred primordial abundance 
adopted here, Yp = 0.241 ± 0.004, by nearly 2a. However, as discussed in §3.3, the 
uncertainties in Yp arc likely dominated by systematic, not statistical errors, so it is 
difficult to know if this tension between D (and '^He and the CBR) and ''He is cause 
for serious concern. In contrast to the other light nuclides, the BBN abundance of 
^He is insensitive to the nuclear reaction rates and, hence, to the nucleon density 
at BBN. Yp is largely set by the neutron to proton ratio at BBN, so that ^He is 
an excellent probe of the weak interactions and of the early Universe expansion 
rate. Perhaps the ^He challenge to SBBN is a signal of new physics. The SBBN 
predicted abundance of ''Li ([Li]p = 2.65;to og) nearly a factor of two higher than 
the observationally determined value^' ([Li]p — 2.37 ± 0.05). While there is some 
spread in the lithium abundances inferred from the data^^, the largest cause for 
concern is that Li is observed in the oldest stars in the Galaxy, which have had 
ample time to modify their original surfaces abundances. For ''Li, it appears likely 
that astrophysical uncertainties dominate at present. 

Thus, while there appears to be qualitative confirmation of the standard mod- 
els of cosmology and particle physics extrapolated back to the first seconds of the 
evolution of the Universe, precision should not be confused with accuracy. The 
accuracy of the presently-inferred primordial abundances of D, "^He, ^He, and '^Li 
remains in question and it would not be at all surprising if one or more of them 
changed by more than the presently-quoted errors. After all, there are only 5 (6) 
lines of sight where deuterium is observed in high-redshift, relatively unprocessed 
(low metallicity) material; ^He is only observed in the chemically processed in- 
terstellar medium of the Galaxy and the lack of variation of its abundance with 
metallicity or with position in the Galaxy STiggcst a very delicate balance between 
post-BBN production, destruction and survival; systematic errors and corrections 
to the ^He abundance inferred from observations of low metallicity, extragalac- 
tic Hii regions are likely larger, maybe much larger, than the current statistical 
uncertainties; lithium is derived from observations of very old, very low metallicity 
stars (good!) in our Galaxy (bad?) and the corrections for stellar atmosphere models 
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and, especially, for main sequence mixing-induced depletion and destruction remain 
large and uncertain. Much interesting, important work remains for observational 
and theoretical astronomers. 

The current standard models receive strong support from these messengers from 
the early universe, confirming in broad outline our understanding of the evolution of 
the Universe and the particles in it, from the first seconds to the present. Any models 
of new physics must consider this success and avoid introducing new conflicts. Much 
interesting, important work remains for cosmologists and high energy theorists. 
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